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I. INTRODUCTION 



The mathematical problem of spanning trees on a connected graph can be considered 
as a model of statistical mechanics and, as such, it is the first non-trivial exactly solved 
multidimensional problem thanks to the famous Kirchhoff's theorem l|. In modern classi- 
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fication, the model belongs to the class of free-fermion models [2| which admit determinant 
solutions. The spanning trees are associated with a variety of models, such as the Abelian 



sandpile [3|, Hamiltonian walks on the Manhattan lattice and dense polymer models 
The enumeration of spanning trees on the two-dimensional square lattice is equivalent to the 

Q n 

close packed dimer problem solved by Kasteleyn [7| and Temperley and Fisher [8|. In the 
scaling limit, correlation properties of the s pan ning trees can be described by the conformal 
field theory with central charge c = —2 9|- 12l |. 

The simplest generalization of spanning trees are the spanning webs, spanning subgraphs 



of a connected graph containing cycles together with tree branches attached to them |13 |. 
There are two sources for the appearance of topologically different classes of cycles in models 
associated with spanning trees. First, periodic boundary conditions in at least one spatial 
dimension generate cycles. Such cycles appear in the exact solutions of the dimer problem 
on lattices wrapped on a cylinder or torus [7| and they are non-contractible to a point in the 



embedding surface. Second, lattice defects like monomers in a dense dimer packing give rise 

nn 

to a different kind of contractible cycles |14l . I15| . In general, the spanning web model does 



not belong to the free fermion class and, moreover, it is not exactly solvable. However, for 
particular geometries of the cycles and appropriate statistical weights of the configurations 
it retains the free-fermion properties. 

In this paper we calculate the partition function of a spanning web model on a finite 
cylinder by considering the number of cycles winding around the cylinder as a parameter. 
Our aim is to evaluate the leading finite-size corrections to the free energy in the limit of 
large perimeter of the cylinder. In the absence of cycles the finite-size effects of the spanning 
tree model confirm predictions of the logarithmic conformal field theory [l^. We show that 
the presence of cycles changes the Casimir effect in accordance with conformal weights 
which appear in the Kac table l^] • An example of exactly solvable logarithmic models with 
conformal boundary conditions has been given recently by Pearce and Rasmussen . They 
considered critical dense polymers with certain types of defects on a strip and reproduced 
the conformal weights in the first column of the extended Kac table. Their results were 
obtained by means of a functional equation for commuting transfer matrices, formulated in 
terms of the planar Temperley-Lieb algebra. 

Despite the similarity between dense polymers at the free fermion point and spanning 
trees on an auxiliary sublattice, the classification of conformal weights in these two models 
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is quite different. The entries of the Kac table for the model of dense polymers are labeled 
by the number of "defect" lines 1^ under fixed boundary conditions at both sides of the 
strip. We will show that the cycles in the spanning webs play the role of pairs of defect lines 
in the model of dense polymers. However, in our case the boundary conditions are different 
for the odd and even entries in the first column of the Kac table. Some complication of the 
boundary conditions is the fee to be paid for simplicity of derivation of the partition function 
of the spanning web model. Analytical calculations in Sections II and III are reduced to 
the standard determinant expressions for the free fermion model with subsequent analysis 
by use of the Euler-Maclaurin formula. In Section IV we show that the partition function 
calculated for a finite lattice with different boundary conditions coincides with the characters 
of coinvariants calculated in different modules for the algebra of symplectic fermions. This 
allows us to identify open and closed boundary conditions with modules generated by integer 
and half-integer modes of fermions respectively. 

II. THE SPANNING WEBS MODEL 

We consider the labeled graph Q = {V, E) with vertex set V and set of bonds E. The 
vertices are sites of the square lattice Sx,y, 1 < x < M, 1 < y < N from which we obtain 
a graph on a cylinder by identifying Sx,y and Sx+M,y for all x, y. The graph Q represents 
a finite square lattice embedded in a cylinder of height and perimeter M, with closed 
boundary conditions at the top and bottom edges. The term 'closed' means the absence of 
bonds connecting vertices of V with an exterior of Q. We shall consider also the case of open 
boundary conditions at the vertices B G V belonging to one of the edges {s^,!, 1 < x < M} 
and {sx,N, 1 < a; < M} of the cylinder, or to both of them. These cases correspond to a 
graph Q' = [V, E') with vertex sei V' = V VJ g containing an additional vertex, the root 
and the set of bonds E' = EU{{j, g) '■ j G B} enlarged with the bonds connecting the vertices 
of B with the root g. For convenience of notation we label the boundary conditions by the 
superscript {fi,i^): /i = (/i = 1) denotes closed (open) top boundary and u = {u = 1) 
closed (open) bottom boundary, respectively. We find it convenient to construct the desired 
spanning web configurations on the above graphs by using the arrow representation, see e.g. 



191]. Accordingly, to each vertex i G V we attach an arrow directed along one of the bonds 



incident to it. Each arrow defines a directed bond {i i') and each configuration 
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of arrows A on Q defines a spanning directed graph (digrapfi) QsdiA.) with set of bonds 
Esd{A) = {{i ^ i') : G V} depending on A. Similarly, the arrow configurations on Q' 
define a spanning digraph Gg^iA) with set of bonds E'^^{A) = {{i ^ i') : i & V, i' & V U g}. 
Note that no arrow is attached to vertex g, thus it has out-degree zero. A cycle of length k 
is a sequence of directed bonds (zi, 12), (^2, "^3), (^3, u), • • • , {h, h) where all ij, 1 < j < k are 
distinct. If both (i i') and {i' i) belong to the spanning web we say that it contains a 
cycle of length 2. Our aim is to study sets of spanning digraphs with no other cycles than 
those which wrap the cylinder. The relevant configurations will be enumerated with the aid 
of a generating function defined as the determinant of an appropriately constructed weight 
matrix. 



A. Cylinder with closed boundaries 

We begin with the examination of the determinant expansion of the usual Laplace matrix 
A for the graph Q. Let the vertices of the set V be labeled in arbitrarily order from 1 to 
n = \V\ = MN. Then A has the elements G {1, . . . , n}) 



Zi, if i = j, 

-1, if i and J are adjacent, (1) 
0, otherwise. 



where Zi is the order of vertex i. Since the matrix A has a zero eigenvalue, its determinant 
vanishes. On the other hand, the Leibniz formula expresses the determinant of A as a sum 
over all permutations a of the set {1, 2, . . . , n}: 

det A = ^ sgn(cr) Ai,^(i) A2,<x{2) • • • \i,a{n) = 0, (2) 

where Sn is the symmetric group, sgn(cr) = ±1 is the signature of the permutation a. The 
identity permutation cr = 0";^ in ([2]) yields the term 2:1 2:2 ■ ■ ■ Zn equal to the number of all 
possible arrow configurations on Q. 

In general, each permutation a G Sn can be factored into a product (composition) of 
disjoint cyclic permutations, say, o" = Ci o C2 ■ ■ ■ o c^. This representation partitions the set 
of vertices V into non-empty disjoint subsets - the orbits 0j of the corresponding cycles 
Cj, i = l,...,k. More precisely, if 0j = ^4,2, • • • ^^i,/,} C ^ is the orbit of Cj, then 
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U^=i0i = V and Yl^=i^i — '^^ where U is the cardinahty of the orbit 0j, equivalently, the 
length of the cycle Cj. The orbits consisting of just one element, if any, constitute the set 
Sfp{a) of fixed points of the permutation: 5'/p(cr) = {v = o-{v),v G V}. In the case of the 
identity permutation aid G S'„ all orbits consist of exactly one element, 0j(crid) = {vi} C V, 
i = 1, . . . ,n, and Sfp{(Jid) = V. A cycle q of length \ci\ = k > 2 will be called a proper cycle. 
A proper cycle of length 2 corresponds to two oppositely directed edges which connect a 
pair of adjacent vertices: (fj^i — > ^^4,2); (^i,2 Note that the vertices of an orbit 0j of 

cardinality = |Oj(a") | > 3 are connected by a closed path on Q which can be traversed in two 
opposite directions: if q is the cycle defined by Vi^i cr(f.t,i) = Vi^2, ^ ■ ■ ■ ^ cr(t;j = 1, 
then the reverse cycle c- can be represented as Vi^. — > cr(fj = Vi^.-i, ^ . . . ^ = '^iM- 

Now we take into account that the proper cycles on Q are of even length only, hence, the 
signature of every permutation in the expansion of the determinant depends on the number 
of proper cycles in its factorization, i.e., if a = Ci o C2 ■ ■ ■ o Cp, where |cj| > 2, i = 1, . . . ,p, 
then sgn((7) = (—1)^. Thus, the terms in Eq. ([2]) can be rearranged according to the number 
p of disjoint proper cycles as follows: 

n [n/2] p 

i=l p=l cr=cio---oCp i=l j^Sfp{a) 

Here is the /c-fold composition of the cyclic permutation q of even length Vi E 0i(cr), 
so that c'^~^{vi) 7^ c^(fi) and c''(wj) = Vi. Note that all non- vanishing off-diagonal elements 
are equal to —1. 

The above expansion reveals the following features: (i) As expected, all spanning digraphs 
on Q have at least one proper cycle; (ii) Each term on the right-hand side with Sfp{a) ^ 
represents a set of HjeSj distinct spanning digraphs which have in common the specified 
cycles ci,...,Cp, and differ in the oriented edges outgoing from the vertices j G Sfp{a). 
These oriented edges may form cycles on their own which do not enter into the list Ci, . . . , c^; 
(iii) Since the sets Uf^^0i and Sfp{ci, . . . ,Cp) are disjoint, the proper cycles formed by the 
oriented edges incident to the fixed points of a given permutation a = Ci o C2 o ■ ■ ■ o Cp should 
enter into the enlarged list of cycles Ci,C2, . . . ,Cp, . . . ,Cp', p' > p, corresponding to the cycle 
decomposition of another permutation a'. 

For example, consider the determinant of the Laplacian matrix of a cylinder of height 
3 and perimeter 4. In the case of closed boundary conditions, the corresponding Leibniz 
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expansion will contain the term 

(-l)3(Ai,5A5,i)A2,2(A3,7A7,8A8,4A4,3)A6,6(A9,i2Ai2,iiAn,10Ai0,9), (4) 

which represents, up to the sign, 22^6 = 12 spanning digraphs on Q with 3 specified cycles 
and all possible oriented bonds outgoing from the vertices 2 and 6, as illustrated in Fig. [1] 
(a). One of these 12 spanning digraphs will contain the additional cycle 2 ^ 6 ^ 2. In the 
case of open boundary at the bottom and closed at the top, the term 

(-l)2(Ai,5A5,l)A2,2A3,3A4,4(A6,ioAio,9A9,i2Ai2,llAii,7A7,6)A8,8 (5) 

represents 22^3^42^8 = 256 spanning digraphs on Q' with 2 specified cycles and all possible 
oriented bonds outgoing from the vertices 2, 3, 4 and 8, as illustrated in Fig. [1] (b). The 
latter bonds will generate three digraphs with one additional cycle of length 2: 2 — 3 — 2, 
or 3 ^ 4 ^ 3, or 4 ^ 8 ^ 4, and one digraph with two additional cycles, 2 ^ 3 ^ 2 and 
4^8-^4. 



As noticed first in lOj, the expansion ([3]) parallels in form the inclusion- exclusion principle 
in combinatorial mathematics. Indeed, let ci,C2, ■ ■ ■ ,Cm be the list of all possible proper 
cycles on Q, labeled in an arbitrary order. Define Ai, i = 1,2, ... ,771 as the set of all 
spanning digraphs on Q containing the particular cycle Cj. Then, expansion ([3]) can be 
written in the form of the inclusion-exclusion principle: 

m 

\uZ,A,\ = J2\A\- E \AnA^nAk\---+{-ir+'\A,n---nArn\, 

i=l ^<i<j^rn l<i<j<k<m 

(6) 

which holds for any finite sets Ai, A2, . . . , A^, where \A\ is the cardinality of the set A. 

Now we are in the position to define a matrix D^^'^\ associated with the graph Q, such 
that det-D(°'°^ be the generating function of all spanning digraphs on Q which have no 
contractible cycles. The elements Df-'^\i,j = 1, . . . ,n = NM of are explicitly given 

as: 

zu if i = h 



r,(o,o) 



— 1, if i,j are vertical neighbors, 

—a, if z is left neighbor of j, (7) 

-a~^, if i is right neighbor of j, 

0, otherwise. 



(a) (b) 

FIG. 1: Possible spanning digraphs on a cylinder, generated by a single term in the determinant 
expansion of the corresponding Laplace matrix (see text). Part (a) corresponds to closed boundary 
conditions, and part (b) to open at the bottom and closed at the top ones. 

Here a = ci;^/*^e~*'^/*^, the terms "left" and "right" fix the opposite directions of the 
horizontal edges. Note that all closed paths which do not wrap the cylinder contain an 
equal number of horizontal edges with either orientations, hence, their weight in det Z}(°''^) 
remains the same as in det A. Therefore, all the configurations which contain such closed 
paths (contractible cycles) cancel out in the expansion of det D^'^'^\ On the other hand, cycles 
generated by off-diagonal elements that wrap the cylinder change their sign, because they 
contain horizontal edges oriented in one direction exceeding by M the number of edges in the 
opposite direction. This amounts to the total factor of a^'^ = —u, or a~^' = —u"^ depending 
on the orientation. Therefore, each non-contractible cycle with a given orientation is counted 
twice, however, with different weight - once it enters into the determinant expansion with unit 
weight, being generated by diagonal elements of the matrix D^^'^\ and second time it enters 
with a factor u or (depending on the orientation) as generated by off-diagonal elements 
of that matrix. Thus, the total number of non-contractible cycles, irrespective of their origin 
and orientation, is given by the coefficient in front of the corresponding power of 00+00^^+2 = 
^ in the series expansion of the partition function. The power of ^ is the "good quantum 
number" which is a well defined and conserved quantity under the action of the transfer 
matrix along the cylinder. However, as we shall see below, more convenient expansions of the 
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partition function, which can be directly compared with characters of the Virasoro modules, 
are given by the power series in uj itself, or in terms of combinations like [uj + uj^^Y and 
Sfc=o'^*~^'^- general, besides the non-contractible cycles, the spanning digraph contains 
tree subgraphs connected to them. All branches of the trees can be generated only by the 
diagonal elements of D^^''^^ and, hence, carry unit weight. 



B. Cylinder with one closed and one open boundary 



The Laplacian matrix A* for the graph Q' corresponding to (/i, z/) = (0,1) boundary 
conditions is a (n+ 1) x (n+ 1) matrix of the same form ([1]) as far as the notions of degree of 
a vertex and adjacency are understood in the context of Q' . However, to make the similarities 
and dissimilarities with the former case apparent, we retain the notation Zi for the degree 
of vertex = 1, . . . , n, with respect to and explicitly introduce the label g for the root 
with degree M in Q' , see Fig. ([2]) b. Thus, for the matrix elements of A* we have 





if 


i = jeV\B, 


Zi + 1, 


if 


i=j eB, 


M, 


if 


i=J=9, 


-1, 


if 


i,j G V are adjacent in 


-1, 


if 


i e B,j = g or i = g,j e B 


0, 




otherwise. 



(8) 



Here B is the set of bottom boundary vertices adjacent to the root g in Q'. Now we make 
use of the fact that by Kirchhoff 's theorem the number Ngt of spanning trees on Q' is equal 
to any cofactor of A* and chose for convenience the cofactor Cgg of the element A*^. Then 



N. 



St 



Cgg = det A', where A' is the n x n matrix with elements = 1, . . . , n): 



A', 



2;^, if i=jeV\ B, 

Zi + 1, ii i = j ^ B, 

— 1, if i,jEV are adjacent, 
0, otherwise. 



(9) 



By comparing the above expression with ([T]) one sees that the only difference is in the 
diagonal elements: the order of the vertices belonging to the open boundary has been in- 
creased by one. Therefore, the same arguments as in the previous section lead us to the 
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c) ()-^ — e-^ — o o Q c*-^ 



tB-^ — — e-^ — o- 



C)-^ — e-^ — e-^ — e-^ y 



C) o — ^-o o — ^-G 



6-^ — e-^ — e-^ — e- 



o — ^—o C) — ^-G 



y 



(a) 




D 



(0,1) 



FIG. 2: Spanning digraphs on a cylinder of height = 6 and perimeter M = 10 under: (a) Closed 
boundary conditions. The two closed paths x — x and y — y represent non-contractible cycles, (b) 
Open boundary conditions at the bottom and closed ones at the top. The three bonds incident 
with the root give rise to boundary trees. 

matrix D^'^'^^ with elements 

Zi, ifi=jeV\B, 
Zi + 1, if i = j e B, 
—1, if i, j are vertical neighbors, 

(10) 

—a, if i is left neighbor of j, 
—a~^, if i is right neighbor of j, 
0, otherwise. 

where, as before, a = cu^/^^e"""^/*^. By construction, in the expansion of det D'-'^'^'' all the 
arrow configurations with contractible cycles cancel out. Thus det Z}(°'^) is the generating 
function of all the spanning digraphs which are either spanning forests of trees rooted at the 
open bottom boundary, or contain non-contractible cycles wrapping the cylinder and tree 
subgraphs rooted at these cycles or at the open boundary. As in the case of closed boundaries, 
each non-contractible cycle with a given orientation is counted twice with different weights 
- once with unit weight, being generated by diagonal elements of the matrix D^^'^\ and 
second time with a factor u or uj~^ (depending on the orientation) as generated by off- 
diagonal elements of that matrix. 

The case of open-open boundary conditions, (/i, z/) = (1, 1), can be considered analogously 
to the previous two cases and we do not describe it separately. 

The partition function calculated as the determinant of the Laplacian, Z^j^^ = detD(^'^\ 
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can be split into a product of two parts, 

4^^^ = e^-4"^Hg,^) (11) 



where e nm is the nonuniversal part of the partition function, including bulk and boundary 
free energy F^^j\ and Z^''^\q,uj) is the universal part which is a polynomial in the aspect 
ratio parameter q = and a Laurent polynomial in to. The universal part of the partition 
function can be decomposed into different combinations of u, for example, 

= f2^[T^^N]{q){u; + uj~'r, (12) 

s=0 
N 

= C^{r\N]i<lP' , (13) 

s=~N 
N s 

Zi!;-\q,u) = Y.Cb'^l^]i^)Y.^""'^ (14) 

valid for all boundary conditions under consideration, /i, z/ = 0, 1. Note that the second and 
third polynomials can be obtained from the first one by using the relations: 

d>:f^[N] = d^l1[N]= j2 ('\^'')ci::%pWl s = 0,...,N, (15) 

p=Q \ P / 

d>:f [N] = d^f [N] -d{;%[N], . = 0, . . . , AT - 2, (16) 
d^f\N] = d^f[N], s = N-l,N. (17) 

Each of the polynomials C,- [A^] (g) , i = 1,2,3, has a well defined statistical meaning. 
For example, C[^g'^\N]{q) is the universal factor in the partition function which is propor- 
tional to the number of configurations with s non-contractible cycles of either orientation, 
generated by off-diagonal elements of the matrix D^^''^^; these configurations may have also 
any allowed (by the size of the cylinder) number of non-contractible cycles of both orien- 
tations, generated by diagonal elements of that matrix. On the other hand, the universal 
factor C^^^^ [N] (q) is proportional to the number of configurations with fixed to s difference 
in the numbers of non-contractible cycles with positive and negative orientation, generated 
by off-diagonal elements of the matrix Z}(^''^); the total number of such cycles, as well as the 
number of non-contractible cycles of either orientation, generated by diagonal elements of 
the same matrix, may take any allowed values. 



11 



In the next section, we shall evaluate the asymptotic form of the coefficients Cf'^^ [N\, 
i = 2, 3, for M ^ oo, ^ oo and disclose their relation to the finitized characters of the 
logarithmic Virasoro modules. 



III. CALCULATION OF THE PARTITION FUNCTION 

The matrices of edge weights D^^''^^ can be written in the form of a sum of direct products 
of simple N x N and M x M matrices: 

D^'^''^) = {2En - Q^i;'''^) ® Em + En® {2Em + ad + a~'Gl), (18) 

where En (Em) is the unit N x N {M x M) matrix, Q^ji^'"^ = {q^'^'''^} is a tridiagonal matrix 
with elements dependant on the boundary conditions [i, j = 1, . . . , N): 

(iff^ = + + + ^i,NSj,N, (19) 

^S'^ = = + + ^^..-1' (20) 

gfy) = + (21) 

Note that 2En - Q^'"^ differs only by sign from the one-dimensional discrete Laplacian 
on a chain of sites with the analogues of Neumann- Neumann (/i = z/ = 0), Neumann- 
Dirichlet (/i = 0, = 1 or /i = 1, = 0) and Dirichlet-Dirichlet (/i = 1, z/ = 1) boundary 
conditions. Hence, the eigenvalues of 2En — Q^-''^^ ^-re X^^'^\p) = 2 — 2cos^, X^^'^\p) = 
ASi'°)(p) = 2 - 2 cos and \^i''\p) = 2 - 2 cos where p = 0,...,N -1. The 

M X M matrix Gri = {gm,n} in ( fT8l) has the elements gm,n = ^m-i,n + ^m,A/^n,i, with m,n = 
1, . . . , M, and is the transposed of Gm- Their eigenvalues are ^luik) = exp{i27ik/M) and 
Uuik) = exp{—i27ik/M), respectively, where = 0, . . . , M — 1. Thus, for the corresponding 
partition functions z\t^^^ = det D^j^'^j we obtain: 

N-l M-1 

^NM = n n [^n"\p) + 2 - ^VMgi.(2fc+l)/M _ ^-l/A/^-i.(2/c+l)/Mj _ (22) 
p=0 fc=0 

By analytic continuation of the identity 

M-l 

+ _ 2 cos ( : 

M 



n 

fe=0 



+ Q-^ - 2 COS ( ^ + a 



g2M^g-2M _2cos(jvfa), (23) 
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with 



g ^ q!(^''^) (p) = 1 + sin^ 4'^) ip) + sin 4'^^) (p) , 

7r(2p + l) (1,1) 



(24) 



2(2iV+ 1 

from real a to complex a = [tt — ilnu;)/M, we obtain 



n{p + 1) 
2(Ar + l) 



p=0 



(UJ + UJ 



-')[Q^^''\p)r'' + [Q^^'^^Hp)]-'''} ■ (25) 



The above exact partition function is a polynomial in x = + uj~^ of degree A^, see (fT2l) . 

When M,N oo, so that M/N = 0(1), a standard finite-size analysis of the contribu- 
tion to the free energy from the u;-independent factor in the partition function yields the 
asymptotic expansion 

N-l 



p=0 



AG 



-MN + M 



TT 



2G 



TT 



ifi + u)- ln(l + V2) 



ttM 



1 
12 



(26) 



where G is Catalan's constant. 




FIG. 3: Correspondence between spanning webs and dense polymers, (a) The cycle x — x is 
surrounded by two polymers, 1 — 1 and 2 — 2, and the cycle y — y hy another two polymers, 3 — 3 
and 4 — 4, which are considered as defect lines in the classification of 



13, 



Iq . (b) Besides the four 



defect polymer lines, there is a fifth line 5 — 5 separating the boundary trees from the rest of the 
lattice. 
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The leading-order asymptotic form of the cu-dependent factor in the partition function 
follows from the approximation Q^^''^^ (p) ^ 1 + ^^''^^ (p) . Thus we obtain 

N-l N-l 



p=0 



p=0 



(27) 



where x = uj + u ^. 

For and M large and N/ M fixed, we rewrite (125|) in the form (fTTl) with the free energy 



Fi^^;) = —MN + M 



2G 



— (/i + I/) - ln(l + 72) 



TT 



(28) 



and universal part of the partition function 

N-l 



3=0 



(29) 



In the next section, we show that fl29l) converges as ^ oo to the characters of symplectic 



fermions 



2l| or, equivalently, to the characters of the doublet algebra A{2) 20|] in the p = 2 



logarithmic model. As is well known, the characters contain complete information about 
conformal dimensions and their multiplicities in the model. 

Before considering the symplectic fermions, we calculate the coefficients C'^^^^ [N] (q) (fT6|) 
and C^^g'^\N]{q) f|T7|) explicitly. The universal part of the partition function fl29|) can be 
rewritten by using the Newton g-binomial formula 

N-l N 

S(3-1) / iV \ 



p=0 s=0 



(30) 



with g-binomial coefficients 



N\ 



N-s+l-^ 



[i-q)---{i-r 



when < s < A^, 



0, when s < or s > N, 



(31) 



in the form (fT3!) with 



N-s 
k=0 



+k{s-l+fi+u) 



N 
s + k 



, < s < A^, (32) 



andC(^/)[Ar](g)=CH[^^](9). 
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For closed-closed and open-closed boundary conditions the summation can be performed 
explicitly and the above expression simplifies to 



2N 
N + S-1 



2N 

N + s 



and 



^^Z/ 2N 
N + s 



(33) 



(34) 



respectively. 

Now we give explicit expressions for C^^/)[iV](g), which are shown in the next section to 
converge to the Virasoro characters as ^ oo. Using f|T6l) . we obtain for < s < — 2: 



- N-s 



1 (fj.-vr I s{s-i+fj.+v) 

ql2 8 ^ 2 



k=0 



N 
s + k 



N-s-l 

E 

k=l 



k'^+k(s-l+tJ.+u) 



N 



N 



(35) 



We note that 

.(0,0) 



(a + l)(s + 2) / 2N — 1 

N + s + 1 



s(s^i) f 2N -I 
q 2 



N + s + l 



s(s + l) /2N - 1 

^ ' \ N + S 



2N-1 
N + s + 2 



(2Ar-l) , 



where Xs^^ (?) ^ire the finitized characters of the logarithmic Virasoro modules, given by Eq. 
5.5) in Ref. 18| for s odd. It should be mentioned that fl36|l resembles Eq. (7.25) in Ref. 



181] for the finitized character of the logarithmic Virasoro module. 
Similarly, 



2N 

N + s 



2s+2 



2N 
N + s + 2 



xZt\<l) (37) 



coincides with the finitized Virasoro characters X2s+2^\'i) given by Eq. (5.5) in Ref. [3] 
for s even. 

On the other hand, under the substitution (A^, s) — > (27V -|- 1, 2s -|- 1) in Eq. (5.5) in Ref. 
[is! ] for s odd, we obtain the relationship 



(38) 
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This relation corresponds to the fact that X^2s+i^\q) is the finitized character of the rank-1 
indecomposable Virasoro module with two irreducible subqotients whose finitized characters 
are given by C^3f2i[N]{q) and C^]s^\N]{q). 
Next, Eq. ((291) implies 



471(9,^) = {oo + 2 + u-')zi^''\q,u;). (39) 
Hence, the coefficients C^f\N] and C^^!.^\n — 1] are related by the equality 

CiTl^M = C^'Mn - m + 2C^];'\n - l]{q) + c£i\[iV - (40) 

which repeats the relations between characters of the logarithmic and irreducible Virasoro 
modules. 

Finally, from Eqs. fl33|) and flM|) it follows that the conformal wights for the open-open 
and open-closed boundary conditions are 

Ai°'°) = ^i^, . = 0,1,2,... (41) 



and 



Af^) = ^^-^, . = 0,1,2,... (42) 



These two sequences can be arranged into the first column of the extended Kac table so that 
its odd entries correspond to (HT!) and the even ones to (l42l) 



It is instructive to compare the obtained results with those of Pearce and Rasmussen Il8 l. 

n 

The partition function of dense polymers evaluated in [18| is a function of the number of 
defect lines / = 0,1,2,... in the polymer system, so that the extended Kac label s = / -|- 1 
runs over all entries in the first column of the table of conformal weights. It is easy to 
notice that even cells differ from their odd counterparts by a fixed value —1/8. This value 
can be associated with the conformal weight hmin of the operator with the smallest scaling 
dimension present in the spectrum of the Hamiltonian, which depends on the boundary 



conditions (see for instance 16| where it is shown that h^[^ = for open-open boundary 
conditions and /imin = — 1/8 for open-closed ones). However, for the dense polymers model 
the boundary conditions at both sides of the infinite strip are equal for even and odd values 
of s. A comparison with the model of spanning webs resolves this illusive contradiction. The 
correspondence between our spanning webs model and the model of critical dense polymers is 
shown in Fig. [3l The left-hand figure shows a spanning web for closed boundary conditions. 
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The polymer lines envelop branches and cycles of the spanning web so that each cycle is 
surrounded by two polymer lines. Then, for closed boundary conditions the number of defect 
lines corresponding to polymers surrounding the cycles is always even. The situation for open 
boundary conditions is shown in the right-hand figure. The branches of the spanning web 
going to the root can be separated from the rest of the web by an additional polymer line 
(marked by 5 in Fig. [3]). Thus, the total number of defect lines becomes odd and one obtains 
the set of even entries of the Kac table. 



IV. CONFORMAL FIELD THEORY OF SPANNING WEBS 



The partition function (129|) has interpretation in terms of symplectic fermions 21|. The 
symplectic fermions are fermionic fields 0^{z) with operator product expansion 

e+{z)e^{w) ^\og{z -w). (43) 

These fields admit periodic and antiperiodic boundary conditions for which they decompose 
with integer 9^, n E Z, and half-integer 9^, n G Z + |, modes, respectively. These modes 
satisfy the anticommutation relations 

[On,^^]+=nSn+m,0- (44) 

Let ^(2) denote this infinite dimensional Clifford algebra. Strictly speaking, ^(2) is not an 
algebra or a vertex-operator algebra because multiplication between integer and half integer 
modes is not defined, but .4(2) is very similar to ordinary vertex-operator algebras and many 



standard notions can be defined for it (see discussion on this subject in [22]). 

The algebra ^(2) contains the Virasoro subalgebra generated by the energy-momentum 
tensor 



T{z) =:d9+{z)d9-{z) : (45) 



with central charge c = —2. 



The algebra A(2) has two irreducible modules Xi and X2 (see details in 20|,l2^). Modules 
Xi and X2 are cyclic with cyclic (vacuum) vectors satisfying = for n > and 

|01) satisfying 6*^101) = for n > | respectively. Module Xi is generated by integer modes 
9^ with n < — 1 from the vacuum vector Module X2 is generated by half-integer modes 
9^ with n < —\ from the vacuum vector |01). The algebra ^4(2) has two projective modules 
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Pi and P2 = X2. The module Pi contains 4 irreducible subquotients isomorphic to Xi 22|. 
The module Pi is cyclic with vacuum vector |00) satisfying ^^|00) = for > 1 and is 
generated from |00) by integer modes 9^ with n < 0. 
For a ^(2)-module X, we define the character 

z) = Tixq^^'^u;'' (46) 

where Lq = ^ zT{z)dz and h is the operator calculating the difference between the 
numbers of 6^ and 6~ modes in a state. 
The character of Xi is 

= n°°a o.) EE^^"''^^(i - (47) 

The character of X2 is 

x'^^^M = tt-71^EE-'"'«'^(i - ,-). (48) 

The character of Pi is 

X^'''\q,u;) = {2 + co + co~')x^'''\q,u;). (49) 
The same characters can be written in terms of Virasoro characters in the form 

r 

X^'''\q,u;) = Y,Y.''''~''Xriiq) (50) 

and 

X^'''\q,u;) = J2i2^'-''Xr2{q) (51) 

where 

^ n-.(i-.V '''' 

are Rocha-Caridi characters of the irreducible Virasoro representations with conformal di- 
mensions 

A„ . (53) 

and c = —2. 
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Now we intend to identify (l29i) with characters of some coinvariants calculated in Xi, X2 
and Pi. (By definition, the coinvariant in the A-module X with respect to the subalgebra 
B (ZAis the quotient X/BX). We fix the subalgebra A{2)[N] of ^(2) for iV G N 

{|6'^„, n > N] periodic 
' - ' (54) 

{6 n> N} antiperiodic 

and consider the characters x^^'^^W]{q^^)^ X^^'^^W]{q^^) ^iid x^^'^H^]iQ^^) coinvariants 
with respect to ^(2)[A^] in the modules Xi, X2 and Pi, respectively. These characters 
coincide with (12^ 

X^^^^\N]{q,u) = Z^^^''\q,uj). (55) 

For N ^ 00 one has 

x(^'^)[iV](g,c.)-.x^^'^H?,^)- (56) 

The polynomials di'^\N]{q) and Cf'^\N]{q) converge to the characters of irreducible 
Virasoro modules 

Cffmq) - XsM. Cffmq) - Xs2{q\ (57) 
and Cff^\N\{q) converges to the characters of logarithmic Virasoro modules for the 



£7V1(1, 2) model 17| as tends to infinity. Note also that characters (!55|) can be expressed 



(1') 

through the Kostka-like polynomials Ki,^{q,uj) from 



in the following way 



X^'''\NUu) = q^^Ki%{q,u;), x^'^'\N]{q,u;) = J^g^+,(g, ^). (58) 

Now we can identify open and closed boundary conditions of the spanning webs model 
with irreducible modules of the algebra ^(2). In conformal field theory, the boundary 
conditions are in one to one correspondence with irreducible modules of the chiral algebra 
23|. For the partition function Z^^''^^ on a cylinder with boundary conditions fx and u the 



conformal field theory predicts 

ZM{q) = Y^N;^X,{q), (59) 

7 

where N'^^ are structure constants in the fusion between modules labelled by fi and u and 
X-yil) are characters of modules appearing in this fusion. The fusion of A{2) modules is 

Xi®Xi = Xi, Xi®X2 = X2, X2®X2 = Pi. (60) 
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The fusion together with identification (l55l) leads to the correspondence 

open b.c. ^ Xi closed b.c. X2. (61) 



An identification of triplet W algebra or Virasoro boundary conditions [2^ is more subtle 
because of their nonlocal nature. We start the consideration with W boundary conditions. 
Each ^(2)-module decomposes into direct sum of two VT-modules, which are labelled by the 
same symbol with additional superscript ± 

Xi = X+©Xf, X2 = X+(BX^, P^ = P+®P{. (62) 

The characters of W modules are obtained from characters of ^(2) modules by taking odd 
and even parts in uj. For example 

xt{q) = \[x^'^'\q.l)±X^'^'\q.-l)\. (63) 

where Xiio) ^i-re the characters of Xf. We can interpret this in lattice terms as taking only 
configurations with even or odd differences between numbers of non-contractible loops of 
different orientations (see also the paragraph after Eq. f|T7|) ). However, it is not clear how 
to formulate such conditions as local boundary conditions without reference to the bulk. 

In order to establish a connection with boundary conditions corresponding to Virasoro 
representations we note that the algebra ^(2) admits a sl{2) action such that O'^^z) and 
9~{z) are the highest and lowest weight vectors of the doublet. This sl{2) action commutes 
with the Virasoro subalgebra (H5l) and therefore irreducible A{2) modules decompose as 

where vr^ is r dimensional irreducible sl{2) representation and Yrs is the Virasoro irreducible 
representation with the conformal dimension (153|) . The character of vr^ is Yl^j=o w^^"^"'? which 
explains the decompositions ( l50l) . (ISTl) and our definition of polynomials C^^g'^\N]{q) in 
f ll4p . A lattice interpretation of C^^g'^\N]{q) is very cumbersome and appeals to conditions 
on non-contractible loops in the bulk like in the W case. 



Conclusions 



In this paper we have found the exact partition function for a model of spanning webs 
parameterized by the number of non-contractible cycles for the finite cylinder geometry. 
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We have calculated the leading finite-size corrections and identified them with the finitized 
characters for the minimal logarithmic conformal field theory with c = —2. 

The model considered here is similar in many aspects to the model of critical dense poly- 
mers solved by Pearce and Rasmussen by using the planar Temperley-Lieb algebra and 
commuting double-row transfer matrices. There are, however, several features which are 
different in these models. First, the cylinder geometry admits classification of the webs con- 
figurations in terms of numbers of non-contractible cycles which are well defined " quantum 
numbers". Regarding the conservation law in the system of dense polymers, one notices 
that the transfer matrices used in have a block triangular structure. This structure 
reflects the fact that defect lines can be annihilated in pairs and, therefore, the number of 
defects is not conserved. Imbedding the system into the cylinder geometry, which is the case 
of our model, is equivalent to taking the trace of a transfer matrix, hence, the number of 
non-contractible cycles (or defect lines in the case of dense polymers) becomes automatically 
fixed. The simple geometry of the cylinder with closed or open boundary conditions on the 
edges allows an elementary evaluation of the partition function by using an extension of the 
Kirchhoff theorem. 

The second and more important feature of the spanning webs model consists in the 
perfect coincidence of the universal part of its partition function Z^''^\q,uj) for different 
combinations of closed and open boundary conditions, /i = 0, 1 and u = 0, 1, with the 
finitized characters of the symplectic fermions, see (|55l) . This allows us to interpret the 
simplectic fermion model as a conformal field theory of spanning webs on a cylinder. The 
further identification of the triplet W algebra in terms of spanning webs is an interesting 
open problem. 

Another problem for future investigation is the explicit construction of Virasoro repre- 
sentations with finitized characters given by Eqs. (36) - (38). Strictly speaking we have 
not proved the indecomposability of these Virasoro representations. However, the structure 
of characters (36) - (38) gives indication that our model belongs to the universality class 
of the c = —2 logarithmic conformal field theory (LCFT). An additional argument in sup- 
port of the above conjecture provides the one-to-one correspondence between the recurrent 
configurations of the Abelian sandpile model (ASM) and the spanning trees on the square 



lattice 



251 . Under this mapping the height variables of the ASM correspond to nonlocal 



correlations in the spanning trees. Explicit calculations of height correlations in the ASM 
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show logarithmic corrections 26(] in complete agreement with the predictions of the c 
LCFT. 
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